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Patterned surfaces with large effective slip lengths, such as super-hydrophobic surfaces 
containing trapped gas bubbles, have the potential to greatly enhance clcctrokinctic 
phenomena. Existing theories assume either homogeneous flat surfaces or patterned sur- 
faces with thin double layers (compared to the texture correlation length) and thus 
predict simple surface-averaged, isotropic flows (independent of orientation). By ana- 
lyzing electro-osmotic flows over striped slip-stick surfaces with arbitrary double-layer 
thickness, we show that surface anisotropy generally leads to a tensorial electro-osmotic 
mobility and subtle, nonlinear averaging of surface properties. Interestingly, the electro- 
osmotic mobility tensor is not simply related to the hydrodynamic slip tensor, except 
in special cases. Our results imply that significantly enhanced electro-osmotic flows over 
super-hydrophobic surfaces are possible, but only with charged liquid-gas interfaces. 



1. Introduction 

The development of microfluidics has motivated interest in manipulating flows in 
very small channels, which exhibit huge hydrodynamic resistance to pressure-driven 



flow (Stone et al. 2004 Squires & Quake 20051. One avenue for driving flow on such 



scales is to exploit hydrodynamic slip, usually quantified by the slip length b (the dis- 



tance within the solid at which the flow profile extrapolates to zero) ( Vinogradova] 1999 
Lauga et aZ.||2007 Bocquet fc Barrat]|2007 1 . For hydrophobic smooth and homogeneous 
surfaces b can be of the order of tens of nanometers (Vinogradova & Yakubov 2003 



Cottin-Bizonne et al. 2005 Joly et al. 2006 Vinogradova et al. 20091, but not much 



more. Since the efficiency of hydrodynamic slippage is determined by the ratio of b to the 



scale of the channel h ( Vinogradova 1995 ) , it is impossible to benefit of such a nanometric 



slip for pressure-driven microfluidic applications. 

In principle, this limitation does not apply to interfacially-driven flows, such as electro- 
osmosis past a charged surface in response to an applied electric field. The combination 
of these two strategies can yield considerably enhanced electro-osmotic (EO) flow on 
hydrophobic surfaces ( Muller et a?.|1986 Joly et a?.|2004 Ajdari fc Bocquet]|2006 ), even 
for nanometric slip lengths. The reason is that the thickness of the electrical Debye layer 
(EDL), characterized by the Debye screening length Ad = defines an additional 

length scale of the problem, comparable to b. For a small surface charge density q, simple 
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arguments show that the electro-osmotic mobility M e , which relates the effective electro- 
osmotic slip velocity (outside the double layer) to the tangential electric field u s = M e E t , 
is given by ( |Muller et alfll986| |Joly et aL||2004[ ): 



M, = 



V 



bn) = - — (I 



(1.1) 



where e and r\ are the permittivity and viscosity of the solution, respectively, and C, = 
q/ne is the zeta potential across the diffuse (flowing) part of the double layer. The 
factor (1 + bn) associated with hydrodynamic slip can potentially enhance interfacially- 
driven flow in microfluidic devices (Ajdari & Bocquet 20061, electrophoretic mobility 
of particles (Khair & Squires 20091, and electrokinetic energy conversion (streaming 
potential) in nanochannels ( van der Heyden et aZ.|[2006 1. 

For this reason, it is attractive to consider electro-osmotic flows over superhydrophobic 
(SH) surfaces, whose texture on a scale L can significantly amplify hydrodynamic slip 
due to gas entrapment (Vinogradova et al. 1995 Cottin-Bizonne et al. 2003) leading 



to effective b of the order of several microns in pressure driven flows ( Ou & Rothstein 
2005 Joseph et al. 2006 Choi et aT| [2006). Equation (1.1 1 with bn > 1 suggests that a 



massive amplification of EO flow can be achieved over SH surfaces, but the controlled 
generation of such flows is by no means obvious, since both the slip length and the electric 
charge distribution on a SH surface are inhomogeneous and often anisotropic. Despite 
its fundamental and practical significance EO flow over SH surfaces has received little 
attention. Recently, Squires (20081 investigated EO flow past inhomogeneously charged, 
flat slipping surfaces in the case of thick channels (h 3> L) and thin EDL (A^ <C L) and 
predicted negligible flow enhancement in case of an uncharged liquid-gas interface, which 
has been confirmed by molecular dynamics simulations (Huang et al. [2008 ) . However, 
this work cannot be trivially extended to the general case of thick EDL (Ad ^> L), where 



improved efficiency of electrokinetic energy conversion is expected (van der Heyden et al. 



2006). For thick EDL, we might also expect anisotropic EO flows transverse to the applied 
electric field, as in the case of rough, no-slip charged surfaces ( AjdaripOOl I. 

In this paper we provide analytical solutions for electro-osmotic flows over weakly 
charged, textured slipping surfaces. We show that the electro-osmotic mobility is gener- 
ally a tensorial property of the surface, reflecting nonlinear averaging of the slip-length 



and charge profiles, and is not trivially related to the hydrodynamic slip tensor (Bazant 
fc Vinogrado^al[2008| . In §|] we give basic principles, formulate the problem and obtain 



general solutions for longitudinal and transverse textures. Effective slip lengths and exact 
solutions for EO velocity over stick-slip stripes modeling SH surfaces are derived in §[3j 
Implications for the use of SH surfaces to enhance EO flows are discussed in §|4] followed 
by concluding remarks in §[5] 



2. General theory 

2.1. Interfacial mobility tensor 



Ajdari (2001) pointed out that linear electrokinetic phenomena are generally tensorial 



in space and showed that microchanncls with both charge and height variations can 
exhibit transverse electrokinetic effects. Here, we ascribe analogous behavior to a 'thin 
interface', whose thickness A is much smaller than the geometrical scale h, by defining a 
tensorial electro-osmotic mobility via u s — M e E, where E is the electric field and u s is 
the effective fluid slip velocity just outside the interface (relative to the surface velocity). 
This is analogous to the tensorial hydrodynamic mobility M^, defined by u s — M^r 
in terms of the normal traction t = n • cr, or equivalcntly, to the slip- length tensor 
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b = Mh?7 defined by u s = b-y for a Newtonian fluid, where 7 = 1-/77 is the strain 



rate (Bazant & Vinogradova 2008). We combine these effects in a general 'interfacial 



constitutive relation' 



M e 



T 
E 



M 



T 

E 



(2.1) 



which acts as an effective boundary condition on the quasi-neutral bulk fluid, where 
K s is a tensorial surface conductivity and j s is the surface current density (in excess 
of the extrapolated bulk current density, integrated over the interface). Using matched 



asymptotic expansions (Chu & Bazant 20071, Eq. (2.1l can be derived by considering 



a semi-infinite quasi-equilibrium electrolyte and solving for the velocity u s and current 
j s 'at infinity' relative to the interfacial thickness, e.g. A = max{Au,I/} for a periodic 
texture of period L or a (non-fractal) random texture with correlation length L. 

Following Bazant & Vinogradova ( |2008 1, we note some basic physical constraints on 
M.. In most cases, we expect Ad to be symmetric, as assumed in (2.1 1, by analogy with 



Onsager's relations for bulk non-equilibrium thermodynamics (Groot & Mazur 19621. 



Indeed, this hypothesis can be rigorously established for Stokes flows over a broad class 



of patterned surfaces (Kamrin et al. 2009 1. Here, we focus on the electro-osmotic mobility 
M e for patterned slipping surfaces by calculating the anisotropic electro-osmotic flow in 



response to an applied electric field, but according to (2.1 I the same tensor also provides 



the 'streaming surface current' j s = M e x in response to an applied shear stress. 

The interfacial mobility Ad must be positive definite for a passive interface, which 
produces entropy and does not do work on the fluid. In general, the second-rank tensors 
Mfe, M e and K s could be represented by 3 x 3 matrices to allow for normal flux of fluid 
(or charge) into a porous (or conducting) surface, driven by normal electric fields (or 
tensile stresses) , but here we will only consider 2x2 matrices in the coordinates of the 
tangent plane to describe impermeable, insulating surfaces. This simplifies our analysis, 
since any symmetric, positive definite 2x2 matrix is diagonalized by a rotation: 



Ml 








COS 6 

—sin 



sin 

cos 



(2.2) 



Once the orthogonal eigen-directions 9 = 0,7r/2 are identified, the problem reduces to 
computing the two eigenvalues, AfJ and M^~, which attain the maximal and minimal 
directional mobilities, respectively. 

2.2. Weakly charged, nano-scale striped patterns 

To highlight effects of anisotropy, we focus on flat patterned SH surfaces consisting of pe- 
riodic stripes, where the surface charge density q and local (scalar) slip length b vary only 
in one direction. In the case of thin channels (h <C L), striped surfaces provide rigorous 
upper and lower bounds on the effective slip (eigenvalues of M/,) over all possible two- 
phase patterns dFeuillebois et a?.||2009l; for thick channels, sinusoidal stripes also bound 



the effective slip for arbitrary perturbations in surface height and/or slip length Kamrin 
et al. (20091. Striped SH surfaces have also been used for reduction in pressure-driven 
flows (Ou & Rothstcin 2005), with a typical geometry sketched in Fig.[T^i corresponding 
to Cassie's state of a roughly flat liquid surface over gas bubbles trapped in wells. By 
symmetry, the eigen-directions of M e , M/j, and K s for a striped surface correspond to 
longitudinal (9 = 0) and transverse (9 = tt/2) alignment with the applied electric field 



or shear stress, so we need only compute the eigenvalues for these cases using (2.2) 



We consider a semi-infinite electrolyte in the region y > above a flat patterned surface 
of period L at y — subject to an electric field E in the x direction. The electrostatic 




Figure 1. (a) Sketch of SH stripes: 6 = n/2 corresponds to transverse, whereas 9 — to 
longitudinal stripes; (b) Situation in (a) is approximated by a periodic cell of size L, with 
equivalent flow boundary conditions on gas/liquid and solid/liquid interface. 



potential is given by 4>{x,y,z) — —Eqx + ip(x,y,z), where ip, is the perturbation due 
to diffuse charge. For nano-scale patterns (L < 1/im), we can neglect convection (Pe = 
(q)EoL/r]KD <g; 1 for a typical ionic diffusivity D), so that ip(x,y,z) is independent of 
the fluid flow. We also assume weak fields \E a \L <C \ip\ and weakly charged surfaces 
{\ip\ -C kT/ze = 25/ z mV at room temperature) for a, z : z electrolyte, so that ip satisfies 
the Debye-Hiickel equation with a boundary condition of prescribed surface charge, 

where, k — Xd 1 — (2z 2 e 2 n oc /eK T) 1 / 2 , is the inverse screening length. In this limit, we 
can neglect surface conduction (which tends to reduce electro-osmotic flow) compared to 
bulk conduction (Du = \] s \/(jEqL <C 1), so we will only discuss the tensors M e and Mf,. 
For transverse stripes, q = q(x) = q(x + L), we expand q(x) in a Fourier series, 

oo 

q(x) = (q) + ^2(A n sm(X n x) + B n cos(\ n x)), (2.4) 

n=l 



where (q) is the mean surface charge, and solve (2.3) by separation of variables, 



t/j - i>(x, y) - W e -m, + J- (A n sm{\ n x) + B n cos{\ n x))eS^\ (2.5) 

£K i^i £ \A 2 + K 



where A„ = 2mr/L. For longitudinal stripes, q = q{z) = q[z + L), the potential ip — 
4>{y,z) has exactly the same form (2.5 1 with x replaced by z. 

The fluid flow satisfies Stokes' equations with an electrostatic body force, 

ryAu = -eA-0V</> + Vp, V • u = 0. (2.6a, b) 

To describe local hydrodynamic slip, we apply Navier's boundary condition 

du 

u(x,0,z) = b(x,z)-g^(x,Q,z), y • u(x, 0, z) = 0. (2.7) 

Far from the surface, u approaches the effective, electro-osmotic 'slip' velocity, 
lim^^oo u(x, y, z) — uq, and the derivatives of u remain bounded. By definition, the flow 
is two-dimensional in the eigen-directions of the surface. 

For transverse stripes, we have u = (u(x, y),v(x,y),0), u(x,0) — b(x)u y (x,0), and 
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v(x,0) = 0. Assuming constant 77 and e and using (2.3 1, we can write (2.6 1 as, 

en 2 1 
Au(x,y) = ip(x,y)S/(f>(x,y) + -Vp, 

where the pressure can be eliminated by taking the curl of both sides, 



(2.8) 



A V x u 



V 



(E x x W). 



Taking another curl of this equation and using incompressibility we obtain 

_ en 2 E sd 2 ip d 2 ip 
rj V dy 2 dz 2 



(2.9) 



where V ip = k ip from (2.3 1. The general solution of (2.9 1 for u(x,y) has the form, 



u(x, y) = u + (P„sin(A„a;) + <3„cos(A„x))e 



K 2 r] dy 2 ' 



(2.10) 



where P n , Q n are unknown coefficients and the last term involving ip(x,y), can be ob- 
tained from (2.5 1. The slip boundary condition then determines the coefficients, {P n , Q n }, 
and the electro-osmotic slip, uo — M^E . 

For longitudinal stripes, the flow is also two dimensional: u = (0, v(y, z), u(y, z)), 
u(0,z) = b(z)u y (0, z), and v(0,z) — 0, where u is again the tangential velocity. Similar 
steps lead to (2.9 1 for u(y, z), using (2.3 1 with ip = ijj(y,z). The general solution now 
takes the form, 



u(y, z) =u + ^2 (PnSin(\ n z) + Q n cos(A„z))< 



-X n y . 



1>- 



(2.11) 



where {P n , Q n } and Uq = mIe o are determined by the slip boundary condition. 



3. Striped super-hydrophobic surfaces 

3.1. Hydrodynamic mobility tensor 

To illustrate the theory, we consider an idealized, flat, periodic, charged, striped SH 
surface in the Cassie state, sketched in figure [TJi, where the liquid-solid interface has no 
slip (61 = 0) and the liquid-gas interface has perfect slip (62 = 00). Let <j>i and </>2 = S/L 
be the area fractions of the solid and gas phases with <fri + <p2 = 1- Our results apply 
to a single surface in a thick channel (h 3> max{Ao, L}), where effective hydrodynamic 



slip is determined by flow at the scale of roughness (Bocquet & Barrat 20071, but not 
to thin channels (h <C mm{Xi>, L}) where the effective slip scales with the channel 
width ( |Feuillebois et q?.||2009| ). 

Pressure-driven flow past stick-slip stripes has been analyzed and shown to depend on 



the direction of the flow ( 


Lauga & Stone 


2003 


Cottin-Bizonne et al. 


2004 


Sbragaglia 


& Prosperetti 


2007 


). Following |Bazant & Vino 


gradova| ( 


2008 1, the hydrodynamic slip 



tensor must have the form (2.2 1, where the eigenvalues are 



M x = = ln 
h 7] 2irn 



sec: 



7T02 

2 



and Ml = 2M£ 



(3.1) 



Note that the effective slip for parallel stripes M'J is twice that of perpendicular stripes, 
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Figure 2. Eigenvalues of the electro-osmotic slip tensor M e for stick-slip stripes of period L 
and slipping area fraction 5/L — 1/2 as a function of electric double-layer (EDL) thickness ac _1 . 
Numerical solutions are compared with exact results in the text for the limits kL — > 0, oo for 
(a) an uncharged slipping interface qi = go, 92 = 0, with thick EDL limit clarified in the inset, 
and (b) a surface of constant charge qi = q2 — go- 



A/jj, analogous to the result from slender body theory that an vertically oriented elon- 
gated body sedimenting due to its own weight falls twice faster then if it were oriented 
horizontally ( jBatchelor] ( |1970[ )). 



3.2. Electro-osmotic mobility tensor 

Using ( |2.2[ ) , we need only calculate the eigenvalues of M e for transverse and longitudinal 
stripes. For our model SH surface with transverse stripes (6 — 0), the region |a;| ^ \8 
has 6 = oo (i.e. u y (x,0) = 0) and q = qi, while the region |<5 < \x\ < \h has 6 = and 



q = q2- Imposing these boundary conditions on the general solution (2.101 yields a dual 
cosine series, 

^ + V(A„Q n + ^^)cos(A n x) = 0, VH<|5, (3.2a) 

' n=l ' 

Mo+ M^ + V(g„+^^)cos(A„x) = 0, V|<5<H<ii, (3.26) 
' n=l ' 

where j n — J A^ + k 2 . (The sine terms vanish due to symmetry.) 

For the general case q\ ^ 02, the dual series can be solved numerically for M^- = uo/Eo 
(see Fig. |2| by truncating the series and taking the inner products with 1 and cos(A n x), 
but exact results are possible in the thin and thick EDL limits. (Below, we also give an 
exact solution for any value of kL in the case gi = g 2 .) In the thin EDL limit, A„/k — ► 0, 
we have j n /K — * 1. Since B n are Fourier cosine coefficients of q(x), the dual series can 
be written as, 



^ A„Q„cos(A„ir) = , V \x\ ^ \5, (3.3a) 

n=l 1 

(n + —J + QnCos(X n x) = 0, V \5 < \x\ < \h. (3.36) 



n=l 
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This dual series can be solved exactly ( Sneddon| [l966 ) to obtain 

r^.thin _ u 91 + 2q 2 K&eff 



In the thick EDL limit, A„/k ^> 1, the dual series (3.2 1 takes the form, 

j n E n B n • 



U 



71=1 

(q)Eo 
r\K 



cos(\ n x) 



Eg 
2r? 



((g) + (ft) V |x| < U 



which can again be solved exactly to obtain the thick-DL electro-osmotic mobility, 



M, 



±,thick _ ^0 
-El 







(g) + 92 

(9) 



(3.4) 

(3.5a) 
(3.56) 

(3.6) 



Rotating the SH surface by 9 — tt/2 for longitudinal stripes, the region \z\ ^ ^5 has 
b = oo and q — qi, while the region |<5 < \z\ |£ has 6 = and q = q 2 . Applying 
boundary conditions to (2.111, we obtain another dual cosine series 



/ E B \ 

]T (A„Q„ + -^^Jcos(A„z) = 0, V \z\ < y, (3.7a) 



w 



(g) J B 0+ g^ , E B, 



cos{\ n z) = 0, V \5 < \z\ s$ l-L, (3.76) 



which can be solved numerically (Fig. 2) or exactly for thin and thick EDL. In the thin 
EDL limit, 7„/k — > 1, the dual series (3.7 1 can be simplified using the q(x) series to 

qi + 2q 2 Kb cS 



obtain again (3.3 1 and thus 



(3.8) 



Therefore, we find that the electro-osmotic mobility te nsor is isotropi c in the thin DL 
limit, Ml hm = Ml hin l, consistent with the examples of |Squires| ( |2008| . In general, M e 
must be isotropic for any flat patterned surface in the thin DL limit, since the effective 
EO slip velocity is equal to the surface averaged EO slip velocity (Ramos et al. 20031, 
and thus always in the direction of Eo. 

In the thick EDL limit, we have A„/7„ — > 1, A„/k 3> 1, and the dual series reduces to, 

E n B„\ . {q)E 



£A n (Q„ + ^)cos(A„z) = -^l, V|*|<§f, 
V ln rj ) rj 



U + + £ (Qn + ^)cos(A„z) = 0, V \6 < \z\ < \L. 

v rjK / v ?77 n / 



E{)B n 



n=l 



which can be solved exactly 



M 



/thick (9) /-. 1 n.T \ 

= (1 + 2rco eff j 



(3.9a) 
(3.96) 

(3.10) 



Since M, 



A^.thick 



|| .thick 



we see that M e becomes anisotropic for thick DL. 
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4. Discussion 

4.1. Uncharged liquid-gas interface 

It is instructive to set the charge to zero on the slipping surface, qi = 0, to describe 
the Cassie state of a SH surface with an uncharged liquid-gas interface. The EO flow 
is then related only to the charge q% = qo on the no-slip liquid-solid interface. First we 
consider the thin EDL limit, where the EO mobility is generally isotropic, as noted above. 



Using (3.8 1 we obtain the simple result of Squires (2008), 



M 



/ hi 



= Mt hm I, where Ml Mn = = - 



I ll i 



qo_ 
r)K 



(g) ; 



(4.1) 



where the EO mobility is the same as for a homogeneous, solid no-slip surface with charge 
density q , regardless of the orientation or area fraction of the slipping stripes. In other 
words, there is no EO flow enhancement due to the slipping regions. As explained by 
Squiresj p008| , the liquid appears to slip on the charged liquid-solid interface by electro- 
osmosis, but without any retarding shear stress or amplifying electro-osmotic flow on the 
uncharged, perfectly slipping liquid-gas interface. 

The flow is anisotropic for any finite EDL thickness, and in this case, there is a simple 
relationship between the electro-osmotic and hydrodynamic mobility tensors in the thick 



(3.11 



EDL limit. Using (3.6) and (3.101 for the EO mobility eigenvalues and comparing with 
we find 



rjK 



M h = -^(I+Kb) 



(4.2) 
In the limit 



which is a natural tensorial generalization of the classical formula (1.1 1. 
kL — > (similar to the thick channel limit), hydrodynamic slip becomes negligible, and 
the EO flow is isotropic and driven by the average surface charge, Ml hlck — > —(q)/r]K. 

From the thick EDL results of this case we see that the EO velocities are even smaller 
than in case of thin EDL. Even if (4.1 1 is used for kL w 1, the EO mobility of the SH 



surface is still smaller than having a constant charge on the surface without any slipping 
regions. This is consistent with the molecular dynamics simulations of |Huang et al.\ 
( 2008 1 for perpendicular stripes for 1 M NaCl solution confined in parallel walls with no 



charge on the liquid-gas interface, which also showed no enhancement of EO flow. Since 
the effective slip length depends on the pattern, however, some flow enhancement may 
be possible for thick EDL, even with an uncharged liquid-gas interface, provided that 
effective slip length is quite large. 

4.2. Charged liquid-gas interface 

The situation is very different if the liquid-gas interface carries some charge. There is 
strong experimental evidence that the air-water interface is negatively charged due pres- 
ence of excess OH~ ions (Takahashi (2005 1 ) supported by the molecular dynamics sim- 
ulations of Zangi & Engberts (20051. Adsorption of ions on the water-air interface has 



also been justified theoretically by Shchckin & Borisov (2005) 



To gain some analytical insight into the possible EO flow enhancement, we consider 
the special case of uniform surface charge q\ = q% = qo- In this case, an exact solution is 
possible for arbitrary EDL thickness and any stripe orientation: 



M e I, 



where M e 



-qo 



l 



Ml 



go 

r)K ' 



(1 + 2b cS n) 



(4.3) 



Clearly, a large EO flow enhancement, similar to that of an isotropic surface ( |1.1[ ), is 
possible with a charged liquid-gas interface (Fig. [2b) . This conclusion holds even if the 
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charge is not homogeneous. For thin EDL kL ~ 10 3 , and <j2 / <Zi ~ 0.1 — 1, t^>2 ~ 0.5, 
the theory predicts EO flow enhancement by factor of 10 — 100. For thick EDL, this 
factor approaches 1, but is offset by the prefactor k _1 in the case of constant charge. 
We conclude that in practical applications of electrokinetic phenomena involving SH 
surfaces, e.g. for microfluidic pumping or energy conversion, electrokinetic enhancement 
is possible with a charged liquid-gas interface, perhaps by an order of magnitude, over a 
homogeneous no-slip surface. 



5. Conclusion 

We have analyzed electro-osmotic flows over patterned surfaces of non-uniform charge 



and local slip length. Unlike the approach of Squires (2008), we have obtained general so- 



lutions for arbitrary EDL thickness, surface charge distribution and slip variation, which 
are qualitatively different from the thin EDL limit. We have shown that the electroki- 
netic response of a patterned slipping surface is generally anisotropic and describable by 
a third rank interfacial mobility tensor (2.1 1, analogous to the tensorial linear response of 



a thin microchannel ( |Ajdari||2001[ ). For stick-slip stripes, we have calculated the electro- 
osmotic mobility tensor and shown that it is not simply related to the hydrodynamic 
mobility tensor (|Bazant fc Vinogradova 2008). In particular, a tensorial generalization 



(4.2 1 of the classical formula (1.1 1 does not hold, except for uncharged slipping regions. 



Our results provide some guidance for the design of SH surfaces for electrokinetic ap- 
plications. For an uncharged liquid-gas interface with thick EDL, our results are closely 
related to those of Lauga & Stone (2003 1 for pressure-driven flows, since EO flow enhance- 



ment is directly set by the hydrodynamic mobility tensor. For an uncharged liquid-gas 
interface with thin EDL, we also confirm the result of iSquires (2008) that there is no 



enhancement of EO flow, but the general response for a SH surface can be quite different. 
Our main conclusion is that a charged liquid-gas interface is required to achieve signifi- 
cant enhancement of EO flow. For a uniformly charged SH surface, we show that the EO 



mobility (4.3) is isotropic and can exhibit large enhancement from hydrodynamic slip, 



possibly by an order of magnitude. 

This research was supported in part by the National Science Foundation, under Con- 
tract DMS-0707641 (MZB). 
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